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This series of papers is devoted to identifying and explaining the properties of strongly correlating 
liquids, i.e., liquids with more than 90% correlation between their virial W and potential energy 
U fluctuations in the NVT ensemble. Paper IV [N. Gnan et ai, J. Chem. Phys. 131, 234504 
(2009)] showed that strongly correlating liquids have "isomorphs", which are curves in the phase 
diagram along which structure, dynamics, and some thermodynamic properties are invariant in 
reduced units. In the present paper, using the fact that reduced-unit radial distribution functions 
are isomorph invariant, we derive an expression for the shapes of isomorphs in the WU phase 
diagram of generalized Lennard-Jones systems of one or more types of particles. The isomorph 
shape depends only on the Lennard-Jones exponents; thus all isomorphs of standard Lennard-Jones 
systems (with exponents 12 and 6) can be scaled onto to a single curve. Two applications are 
given. One is testing the prediction that the solid-liquid coexistence curve follows an isomorph 
by comparing to recent simulations by Ahmed and Sadus [J. Chem. Phys. 131, 174504 (2009)]. 
Excellent agreement is found on the liquid side of the coexistence, whereas the agreement is worse 
on the solid side. A second application is the derivation of an approximate equation of state for 
generalized Lennard-Jones systems by combining the isomorph theory with the Rosenfeld-Tarazona 
expression for the temperature dependence of potential energy on isochores. It is shown that the 
new equation of state agrees well with simulations. 



I. INTRODUCTION 

This is the fifth in a series of papers^— investigating 
the properties of strongly correlating liquids,— i.e., liq- 
uids that have strong correlations between their constant- 
volume equilibrium fluctuations of potential energy, U (t) , 
and virial&i W(t) = —1/3 ^ • V ri t/(ri, rjv) where 
Yi is the position of particle i at time t. As is well known, 
the average virial W gives the configurational contribu- 
tion to the pressure: 



pV = Nk B T + W . 



(1) 



Letting A denote instantaneous deviations from equilib- 
rium mean values, the WU correlation is quantified by 
the correlation coefficient R (with (...) denoting equilib- 
rium average): 



R = 



{AWAU) 



y/((AW)*)((AUy) 



(2) 



Perfect correlation gives R = 1. As a pragmatic def- 
inition we have chosen "strongly correlating liquids" to 
designate liquids that have R > 0.9 in the NVT ensemble 
(constant volume, V, and temperature, T). 

Strongly correlating liquids have simpler physics than 
liquids in general, an observation that has particular sig- 
nificance for the highly viscous phased— Thus it has 
been shown that strongly correlating viscous liquids to a 
good approximation have all eight frequency-dependent 



thermoviscoelastic response function: 



,19-21 



given in terms 



of just one^ 2 - (i.e., are single-parameter liquids in the 
sense of having dynamic Prigogine-Defay ratio-^ 9 . close 
to unity£ i 20 i 22 ). Strongly correlating viscous liquids 



moreover obey density scaling 2 ^— to a good approx- 
imation, i.e., their dimensionless relaxation time f = 
Tp 1 / 3 \JksT /to (where to is the average particle mass) 
depends on density p = N/V and temperature as f — 

Paper I 1 presented computer simulations of 13 differ- 
ent systems, showing that van der Waals type liquids 
are strongly correlating, whereas hydrogen-bonding liq- 
uids like methanol or water are not. Strongly correlat- 
ing liquids include^^ i 5 ' 22 i 29 , for instance, the standard 
Lennard-Jones (LJ) liquid, the Kob-Andersen binary LJ 
(KABLJ) mixture, an asymmetric rigid-bond dumbbell 
model, a seven-site united-atom toluene model, and the 
Lewis- Wahnstrom OTP model. 

Paper 112 analyzed the cause of WU correlations with 
a focus on the LJ potential. The strong correlations were 
related to the well-known fact that an inverse power- 
law (IPL) pair potential, v(r) oc r~ n where r is the 
distance between two particles implies perfect WU 
correlation^, 



AW(t) = jAU(t) 



(3) 



with 7 = n/3. Around the potential energy minimum, 
the LJ potential is well described by an "extended" in- 
verse power-law potential (elPL)^ t>Lj(V) = Ar~ n + B + 
Cr. At constant volume the linear term contributes lit- 
tle to the virial and potential-energy fluctuations: When 
one nearest-neighbor interatomic distance increases, oth- 
ers decrease in such a way that the sum is almost con- 
stant. Thus systems interacting via the LJ potential in- 
herit strong WU correlations from an underlying inverse 
power-law - they have a "hidden scale invariance" 4^ 
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Paper III 3 gave further numerical evidence for the ex- 
planation for strong WU correlations presented in Paper 
II, and theoretical results were given on the statistical 
mechanics and thermodynamics of the hidden scale in- 
variance that characterizes strongly correlating liquids. 
It was also shown that strong virial-potential energy cor- 
relations are present even in out-of-equilibrium situations 
- the hidden scale invariance is a property of the potential 
energy surface, not just of the equilibrium states. 

Paper IV— introduced the concept of "isomorphs" in 
the phase diagram of a strongly correlating liquid. Start- 
ing from a single assumption a number of isomorph in- 
variants were derived. In particular, structure and dy- 
namics were shown to be invariant on isomorphs when 
reduced units are used. 

In the present paper further simulation results sup- 
porting the isomorph predictions are presented using sys- 
tems interacting with the multicomponent generalized LJ 
potential: 

%'(r-«)=i m) (n i ) + i n) (r ij ), (4) 

where (ry ) is an IPL potential acting between the 
two particles i and j: 

«J ) (r«)^S , ( ff W*- (5) 

For systems interacting via a generalized LJ potential, 
a prediction for the shape of the isomorphs in the WU 
phase diagram is derived in section III and demonstrated 
to fit well to simulation results 4^ Interestingly, the iso- 
morph shape depends only on the exponents m and n. 
Thus, e.g., all 12-6 LJ systems have the same isomorphs 
in the WU phase diagram. Finally we briefly present two 
applications of the theory. One tests the Paper IV pre- 
diction that solid-liquid coexistence lines are isomorphs. 
The second application gives an approximate equation 
of state for systems interacting via generalized LJ po- 
tentials; this is arrived at by combining the present the- 
ory with Rosenfeld and Tarazona's expression for the iso- 
choric temperature dependence of the potential energy. 



II. ISOMORPHS 
A. Definition 

We term a microscopic configuration "physically rele- 
vant" if its influence on the thermodynamics and dynam- 
ics of the sustem is not a priori negligible (Paper IV). 
For instance, any configuration with very strong particle 
overlap is physically irrelevant; note, however, that even 
unlikely configurations like transition states are relevant. 

Two state points (1) and (2) with temperatures T\ and 
Ti and densities p\ and pi , respectively, are defined to be 
isomorphic (Paper IV) if they obey the following: Any 
two physically relevant configurations of state points (1) 



and (2), (r^ , ...,rjy ) and (r^ 2 \...,r^), which trivially 
scale into one another, 

„l/3 r (U = p V3 r (2) (i = lt ... iJV ), (6) 

have proportional configurational Boltzmann statistical 
weights: 

e -a( r W,...,rW)/fe B T 1 = Ci2e -U(T^,...,r^)/k B T 2 _ ( ? ) 

Here U(ri, ...,rjy) is the potential energy function and it 
is understood that the constant C\i depends only on the 
state points (1) and (2). 

The property of being isomorphic defines a mathe- 
matical equivalence relation on the set of state points. 
The corresponding equivalence classes are smooth curves 
in the phase diagram termed isomorphs or isomorphic 
curves. 



B. The approximate nature of isomorphs 

Equation ([7]) implies identity of the normalized Boltz- 
mann probabilities for scaled relevant configurations of 
isomorphic state points. As detailed in Paper IV this 
identity implies that several quantities are invariant 
along an isomorph. Examples include the configura- 
tional entropy, the isochoric specific heat, N-particle 
distribution functions in reduced units (in particular, 
the radial distribution function(s)), reduced- unit dynam- 
ics (both for Newtonian and stochastic dynamics), nor- 
malized reduced-time autocorrelation functions, reduced- 
unit transport coefficients, etc. It was further shown in 
Paper IV that a jump from an equilibrium state to an iso- 
morphic state point brings the system instantaneously to 
equilibrium (see also Ref. H3). 

IPL potentials are Euler homogeneous functions, i.e., 
obey C/(Ari, Arjv) = A~ n f7(ri, rjv). Such purely 
repulsive potentials - which do not describe real inter- 
molecular interactions - are the only potentials that have 
100% correlation between virial and potential energy fluc- 
tuations; this reflects the well-known IPL virial-theorem 
identity W — (n/3)U. Likewise, only IPL liquids have 
exact isomorphs; they satisfy Eq. J7| with Cn = 1 for all 
configurations of two state points obeying p[jT\ = p\jTi 
where 7 = n/3. 

Appendix A of Paper IV showed that strongly correlat- 
ing liquids generally have isomorphs to a good approx- 
imation. Consequently, to a good approximation these 
liquids inherit a number of the exact invariants that IPL 
liquids have along their curves given by p 1 /T = Const., 
namely all the IPL invariants that do not rely on the IPL 
identity C12 = 1. Examples include the above mentioned 
thermodynamic, static, and dynamic isomorph invari- 
ants. IPL invariants that are not inherited by strongly 
correlating liquids in general include, e.g., the Helmholtz 
free energy over temperature, the potential energy over 
temperature, and the reduced-unit compressibility. 
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The fact that strongly correlating liquids have iso- 
morphs to a good approximation should be understood 
as follows. Any liquid has curves in its phase dia- 
gram of constant excess entropy curves of constant iso- 
choric specific heat, curves of constant reduced relaxation 
time, curves of constant reduced-unit transport coeffi- 
cients, etc. For a strongly correlating liquid these curves 
are almost identical, and they define its (approximate) 
isomorphs. Points on these isomorphs have approxi- 
mately identical structure and approximately identical 
dynamics as probed, e.g., by normalized reduced-time- 
autocorrelation functions. 

For clarity of presentation below we shall not repeat- 
edly mention that the existence of isomorphs is not exact 
for generalized LJ systems. Accordingly, we shall speak 
of isomorphs as a unique, well-defined concept. In other 
words: all properties of isomorphs derived below in gen- 
eralized LJ systems would be exact if these systems did 
have exact isomorphs, e.g., fulfilled Eq. (JTJ); since this is 
not the case, the properties are in principle - as well as 
practice - approximate. 



C. Generating isomorphs in simulations 
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The structure in reduced units ( 



f = p 1 / 3 



-) is predicted 



to be invariant along an isomorph (Paper IV). The ex- 
cess entropy, S cx = S — Sideaii depends only on structure 
and was shown in Paper IV also to be invariant along 
an isomorph (here and henceforth "excess" refers to the 
quantity in question in excess to that of an ideal gas with 
same density and temperature). In the below simulations 
we generate an isomorph as a set of state points in the 
phase diagram with constant S ex . To change density and 
temperature keeping S cx constant, the following identity 
is used (Paper IV): 



91nT 
dlnV 



dV IT 



rn ( dS cx 
J I 3T 



(8) 



A. (9W\ 
V V dT ) V 



and a Maxwell re- 



where /3° x = 

lation was applied. Utilizing the standard fluctuation 
formulae this leads (Paper IV) to: 



<91nT 
<91np 



(AUAW) 
((AC/) 2 ) • 



(9) 



The procedure applied to generate isomorphs numer- 
ically can be summarized as: 1) An equilibrium NVT 
simulation was performed at a given state point; 2) 7 
was calculated from the fluctuations using Eq. (j9|); 3) 
The density p was changed slightly (of order 1%), and 7 
was used to calculate the corresponding change in tem- 
perature in order to keep S cx constant (Eq. ((9}); 4) A 
simulation at the new state point was performed, and the 
procedure was repeated. 



FIG. 1: a) Radial distribution function (AA particles) for 
three isomorphic state points of the Kob- Andersen binary LJ 
mixture 45 (KABLJ). b) Radial distribution function along an 
isotherm of the KABLJ mixture. A similar picture applies to 
the AB and BB radial distribution functions that are, how- 
ever, slightly less isomorph invariant (Paper IV). 



Figure DJa) demonstrates the isomorph invariance of 
the radial distribution function for the large (A) parti- 
cles of the Kob- Andersen binary LJ mixture^ (KABLJ). 
For comparison, Fig.[T](b) shows the same on an isotherm 
with a similar (but smaller) density change. Clearly, 
structure is to a good approximation invariant on the 
isomorph, whereas this is not the case on the isotherm. 

Taking the logarithm of Eq. ([7]) shows that the po- 
tential energy surface in reduced units (U = U/ksT) 
is the same for two isomorphic state points, except for 
an additive constant. This constant does not influence 
the dynamics, and consequently the dynamics in reduced 
units is the same for the two isomorphic state points. 
Figure [2fa) demonstrates the isomorph invariance of the 
incoherent intermediate scattering function of the A par- 
ticles of the KABLJ mixture. For comparison, Fig. [2](b) 
shows the same property on the T — 0.55 isotherm. Like 
the structure, the dynamics is also invariant to a good 
approximation on the isomorph - both regarding average 
relaxation time and the shape of the relaxation function; 
this is far from the case on the corresponding isotherm. 

The exponent 7 calculated from Eq. (|9]) generally de- 
pends on the state point. Figure [3] shows (as red squares) 
7 along the isomorph of Figs. 1(a) and 2(a). Temperature 
changes by almost a factor of two, but 7 changes less than 
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FIG. 2: a) Incoherent intermediate scattering function of the 
large (A) particles for isomorphic state points (density varia- 
tion 11.6%) of the KABLJ mixture; time is given in reduced 
units: t = tp 1 ^ 3 \fkj^F~pm and the q- vector is kept constant 
in reduced units, b) Same as in a), but along an isotherm 
(density variation 9.5%). 
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FIG. 4: Isochores for four different binary generalized LJ 
mixtures. 12-6 KABLJ is the Kob-Andersen binary 80:20 
mixture 45 using the standard 12-6 LJ potential. 18-6 KABLJ 
is the same system with the generalized LJ potential (m = 18, 
n = 6, Eq. Q) and parameters chosen to place the mini- 
mum of the potential at the same distance and depth as in 
the 12-6 KABLJ system. Likewise, 18-9 KABLJ is using the 
exponents m — 18 and n = 9. Finally 12-6 WBLJ is the 
Wahnstrom 50:50 mixture 46 using the standard 12-6 LJ po- 
tential. Straight lines are linear regression fits. The inset 
shows a zoom on the data for the two 12-6 LJ systems. For 
each system the simulated temperatures ranges from where 



non-exponential relaxation sets in (r a f 
state with two-step relaxation (r a « 10 3 



10 u 
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to a viscous 



5%. Paper IV argued that for a system with good iso- 
morphs 7 depends only on density. This is supported by 
Fig. [3] in which 7 on an isochore is shown to vary only 
~ 1% when temperature is changed by a factor of two 
(black circles). 

It follows from Eq. ([8]) that the slope of isochores in 
the WU phase diagram is given by 7: 
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FIG. 3: 7 calculated from equilibrium fluctuations via Eq. ([9]) 
for the KABLJ system. Red squares: state points belonging 
to the isomorph presented in Figs. [T] and [2] Black circles: 
isochoric state points. Straight lines are linear regression fits. 
Standard errors on slopes are 0.06 (isomorph) and 0.04 (iso- 
chore) . 



f dW 
VdU 



(10) 



The insignificant change of 7 along isochores (black cir- 
cles in Fig. [3]) means that these to a good approximation 
are straight lines in the WU phase diagram. This is il- 
lustrated in Fig. U for four different binary generalized 
LJ mixtures. 



III. SHAPE OF ISOMORPHS IN THE WU 
PHASE DIAGRAM 



What is the shape of an isomorph in the WU phase 
diagram? This section answers this question in generality 
for the multi-component generalized LJ potential. Recall 
that this potential (Eq. (QJ) is the sum of two IPL's. 
Correspondingly, the potential energy and virial can be 
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FIG. 5: Schematic drawing of the WU phase diagram for sys- 
tems interacting via the standard 12-6 LJ potential. Included 
in this figure are also two theoretical isomorphs, plotted using 
Eq. (I28p . The accessible points are those above both dashed 
lines (see text). The part of the isomorphs between the dashed 
lines in the W, U > quadrant are thus not accessible to the 
standard 12-6 LJ (they are accessible to the purely repulsive 
12-6 LJ potential with e\f > 0.) 



expressed as sums of two IPL terms: 



U = U m + U n , U k = (Y. v k ? ( r «) ) • ( n 

\ i>3 I 

For pair interactions the virial is given by£ 



W 



(12) 



where the prime denotes the derivative with respect to 
Tij . From this we get 



m n 

w = -u m + -u n 



(13) 



For any point in the WU phase diagram we can solve 
Eqs. Unj and CE3]) for (U m ,U n ): 



U m = 
Un = 



3W-nU 
m — n 
-3W + mU 



(14) 
(15) 



These equations allow one to determine which regions 
of the WU phase diagram are accessible to the system. 
Adopting the convention m > n, we restrict ourselves to 
the case e\j > in order to have a repulsive core for 
all interactions. This implies U m > for all configura- 
tions, and therefore also for all thermodynamic averages. 
Combining this with Eq. ([13]) implies W > §{/. If all 



interactions have an attractive part (e^ < 0) we have 
U n < 0, which in combinations with Eq. (|15|) implies 
W > ^-jU . As an example, the region of the WU phase 
diagram accessible to the standard 12-6 LJ potential is 
given by W > AJJ and W > 2U (the later inequality be- 
ing relevant when U < 0). See Fig. [5] for an illustration. 



A. Parametric description of isomorphs 

Along an isomorph the structure is invariant in reduced 
units (Fig.QJ. From this it follows that Uk/p k ^ 3 is invari- 
ant for any k along an isomorph (f = p 1 ' 3 /*): 



Uk 

k/3 




(16) 
(17) 



If we let "*" denote a reference point, (W 7 *, L 7 ,), we can 
use Eqs. (0 and <jT5j) to get U*) and find for other 
state points on the same isomorph (where p = p/p*): 



fc/3 



Ut = P k/3 U, 



(18) 



Combining Eq. (JTHJ) with Eqs. (fTTj) and p3|) we obtain a 
parametric description of an isomorph in the WU phase 
diagram: 



U = p m / 3 U* m + p n / 3 U* n , 

W = ^rn /3u:n+ n~n /3u: 



(19) 

(20) 



Once the numbers U^ and U* have been determined from 
a reference state point using Eqs. ([T4"f and (|T5|) . the entire 
isomorph to which this state point belongs is thus traced 
out in the WU phase diagram using Eqs. (|19l) and (|20l) . 
Alternatively, given a collection of isomorphic state point 
generated, e.g., as described in section lie, the predicted 
relation between U, W, and p can be tested by linear 
regression: 



-(ro-n)/3rj* , tj* 
r ^ m ' ^ n 



u 

pn/3 

W_ _ rn ,_ n)/3 n 



(21) 
(22) 



This is done in Fig. |6] for the 12-6 KABLJ isomorph pre- 
sented in Figs. [T]and [5J The potential energies were fitted 
by linear regression by means of Eq. ([2"Tj) . Subsequently, 
the virials were compared to Eq. (|2"2"T) using the parame- 
ters estimated from the potential energies (Eq. ([2"Tj) L i.e., 
without performing a new fit. 

We note here two important consequences of the above. 
1) If the reference state point (J7*,W*) generates the 
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FIG. 8: The "master isomorph": collapsing the four iso- 
morphs from Fig. [7J by scaling with Wo defined as the virial 
on each isomorph when (7 = 0. The points give simulation 
results, the dashed curve the prediction of Eq. 12811 ■ 



FIG. 6: Test of 12-6 KABLJ isomorph (same as in Figs. Q] 
and [2]). Top panel: U fitted to Eq. (pi) using p, = 1 (i.e., 
p = p). Bottom panel: W compared to Eq. (|22p . using the 
parameters and found from U (top panel). 
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FIG. 7: Four different 12-6 LJ isomorphs in the WU phase 
diagram. The points give simulation results, the full curves 
the predictions of Eqs. (|19l) and Q2Q0. Two isomorphs be- 
long to the KABLJ system, the third result from simulations 
of Wahnstrom binary LJ liquid 46 (WBLJ), and the fourth 
(SCLJ) is the single component LJ liquid. 



FIG. 9: Isomorph of the 18-6 KABLJ liquid compared to 
the theoretical prediction (Eqs. (|19|l and (|20j) ). For reference, 
also the isomorph prediction for 12-6 LJ systems is included. 



morph" . 



isomorph (U(p), W(p)), then the reference state point 
{aU^.aW^) generates the isomorph (all(p), a W(p)). This 
means that all isomorphs for a given system have the 
same shape in the UW diagram and can be scaled onto 
one another. 2) The shape of the isomorphs depends 
only on the exponents, m and n. This means that all 12- 
6 LJ systems have the same isomorphs in the WU phase 
diagram - they may be scaled onto a single "master iso- 



Figure [7] shows four isomorphs in the WU phase di- 
agram for three different 12-6 LJ systems4£ Figure [S] 
shows the same isomorphs scaled onto the 12-6 mas- 
ter isomorph, Figure [9] shows an isomorph for the 18-6 
KABLJ system and compares it to a 12-6 LJ isomorph. 
Overall, the agreement between simulation results and 
predicted shapes is quite good. 
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B. Eliminating the density parameter: The 
"master isomorph" equation 

Combining Eq. (JTHJl with Eqs. ([14]) and ([15]) we get 

3W -nU 



~m/3 



Ur, 
U* 



~n/3 = £n 

F rj* 



3W* - ntf, ' 
31F-m<7 
31F* - mil* ' 



(23) 
(24) 



Eliminating density and rearranging, we find an invariant 
for the isomorph: 



(W- f u) m _ {W* -fU* 
(W- §?/)" 



(W* - f £/* 



(25) 



Choosing the reference point with zero potential energy 
(17*, W») = (0, Wq) (this reference point is guaranteed to 
exist if all interactions have attraction), we get: 



(w*y 



(26) 



Wq is a unique number identifying the isomorph. Equa- 
tion ([26]) implies that the isomorph is the solution to 



W m U 



3 W* 



W* 



m/n 



W n U 



3 WX 



(27) 



Here it can be seen directly that - as argued above - for 
a given set of exponents (m, n) all isomorphs have the 
same shape, suggesting the name "master isomorph" ; 
Wq is the parameter that determines the scale of each 
isomorph. For m = 2n (e.g., the standard 12-6 LJ po- 
tential) the solution to this (quadratic) equation is: 



, W 

'w* 



= 1 




7777- (28) 



This has real solution(s) whenever U /Wq > — , where 
equality corresponds to W = 0. Figure [5] plots Eq. 
with n — 6 for two values of Wq . 



IV. APPLICATIONS 



A. Solid-liquid coexistence 




20 40 60 80 

U parameters 



FIG. 10: Test of isomorph predictions for the liquid side of the 
solid-liquid coexistence curve of the single component gener- 
alized LJ potential (n = 6). Left panel: Eq. \21\ , right panel: 
Eq. (|22p : data points are simulation results from Ahmed and 
Sadus,— straight lines are linear regression fits. The lower 
panel plots the fit-parameters from Eqs. (|2ip and (I22p against 
each other, demonstrating that the fits to energies and virials, 
respectively, are consistent. 
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FIG. 11: Same data as in Fig. 1101 plotting potential en- 
ergy versus density (left panel) and virial versus density (right 
panel). Dashed lines: Isomorphic predictions (Eqs. (|19[) and 
(|20p ). using the parameters determined in Fig. [TUJ 



In Paper IV it was argued briefly that Eq. (JTJ) predicts 
that solid-liquid coexistence lines are isomorphs. In more 
detail the argument goes as follows. First one notes that 
an isomorph cannot cross a solid- liquid coexistence line: 
Recall that the isomorph concept relates to equilibrium 
ensemble probabilities. For a pair of isomorphic state 
points all pairs of microscopic configurations related by 
Eq. ([B]) have identical Boltzmann probabilities accord- 
ing to Eq. ([7]). Consequently, the two isomorphic state 



points must also have the same macroscopic phase be- 
havior. In other words, if an isomorph were to cross the 
coexistence line, part of the isomorph would have liquid 
microstates as the most likely and part of it would have 
mixed crystalline-liquid microstates as the most likely. 
This would violate the proportionality of Boltzmann sta- 
tistical weights for scaled states, the property that defines 
an isomorph. 

Given that an isomorph cannot cross the solid-liquid 
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FIG. 12: Solid-liquid phase behavior for the single-component 
12-6 LJ system. Main figure: Filled symbols mark the co- 
existence region in the WU diagram. Blue circles are from 
Ahmed & Sadus^ and red squares are from Mastny & de 
Pablo.H Full curves are the isomorph prediction (Eq. (1281) ) 
with Wq /N = 49.8 and 56.8 respectively. Open circles in- 
dicate results from simulations in the solid and liquid phase 
respectively. Both phases are (very) strongly correlating with 
R > 0.99. Inset shows the solid-liquid phase behavior in a 
p, T-diagram with logarithmic axes. Dashed red lines are the 
isomorph prediction p~" /T — const using 7 found using Eq. Q 
from the simulations of the two phases (open circles) giving 
7 = 5.24 and 5.18 respectively. Full black lines are free fits 
giving 7 = 5.13 and 6.49 respectively. 



coexistence line, the next point is to note that - because 
all state points belong to some isomorph - a curve in- 
finitesimally close to the coexistence line must be an 
isomorph. This argument applies in the W, U phase 
diagram, as well as in the density-temperature or the 
pressure-temperature phase diagrams. In the first two 
cases there is a multitude of isomorphs in the coexis- 
tence region; for instance one particular isomorph will be 
characterized by 53% liquid and 47% crystal, etc. In the 
pressure-temperature phase diagram these "coexistence" 
isomorphs collapse into one. 

The prediction that the melting line, both slightly to 
the liquid side and slightly to the solid side, are iso- 
morphs has a number of consequences and it sheds new 
light on well-known phenomenological melting rules (as 
well as exceptions to these rules for non-strongly cor- 
relating liquids, Paper IV). Thus, for instance the ob- 
servation that the Lindemann melting criterion is pres- 
sure independent for a given liquid follows from the iso- 
morph property of the melting line (whereas the exis- 
tence of a universal Lindemann criterion does not fol- 
low). Indeed, as is easy to show from paper IV the 
isomorph theory implies Ross' "generalized Lindemann 
melting law" from 1969^, according to which the reduced 
configurational partition function of the crystalline phase 
Q* = f dRexp[-P(U(R) - 17(0)] is invariant along the 



melting line (here U is the potential energy, J7(0) is the 
potential energy of the crystal with all atoms located ex- 
actly at their lattice sites, and the tilde denotes reduced 
coordinates. Further well-known melting rules, which the 
isomorph theory explains, include (Paper IV): The re- 
duced viscosity, the reduced surface tension, the reduced 
diffusion constant, and the reduced heat conductivity are 
all invariant in the liquid phase along the melting line. 
Likewise, the reduced-unit static structure factor is in- 
variant along the melting curve. These rules have been 
confirmed for several liquids in both experiment and sim- 
ulation (sec references of Paper V) . 

Like other isomorph predictions following from Eq. (J7J , 
the prediction that solid-liquid coexistence curves are iso- 
morphs is only expected to be exact for systems that 
obey Eq. (J?) exactly, i.e., perfectly correlating systems. 
Indeed, it straightforward to show that in IPL systems 
the phase behavior depends only on p"/ 3 /T, and thus the 
co-existence line is an (exact) isomorph. The interesting 
question now is: how well do strongly correlating sys- 
tems, e.g., generalized LJ systems, follow the prediction? 
Note that we do not expect the liquid-gas co-existence 
lines to be isomorphs (and indeed they are not); the gas- 
phase is not strongly correlating and thus Eq. ([7]) cannot 
be fulfilled in the relevant part of configuration space. 
In contrast, at the solid-liquid coexistence both phases 
involved are strongly correlating for LJ systems (Paper 
I), and we thus in the following test to what degree the 
isomorph prediction is fulfilled. 

Ahmed and Sadus used in a recent paper— a novel nu- 
merical method for determining solid-liquid coexistence. 
They reported results for the generalized LJ potential 
with the repulsive exponent m varying from 12 to 7 
and fixed attractive exponent n = 6 (Ahmed and Sadus 
termed the repulsive exponent n, not to; we keep, how- 
ever, the above notation). Figure [10] shows the state 
points reported by Ahmed and Sadus for the liquid side 
of the solid-liquid coexistence curve compared to the iso- 
morphic predictions as expressed in Eqs. (f2"T| and (|22l) . 
The lower panel shows that the parameters found by fit- 
ting to Eqs. (f2TT) and ([22]) respectively are consistent. 
Figure [TT] shows the same data plotting, respectively, po- 
tential energy (left panel) and virial (right panel) versus 
density. Figure[P2lshows the 12-6 LJ solid-liqud phase be- 
havior in the WU-diagram and in the pT-diagram (inset) , 
including here also data from Mastny and de Pablo^ 

The prediction that solid-liquid coexistence lines are 
isomorphs agrees well on the liquid side with the simu- 
lation results of Ahmed and Sadus 5 ^ as well as with the 
results for 12-6 LJ from Mastny and de Pablc 5 ^. On the 
solid side the agreement is less convincing (see Fig. [T2l : 
in particular the logarithmic derivative in the pT dia- 
gram (slope in the inset) differs by 25%. At present we 
do not have an explanation for this. It might be related 
to inter-particle distances longer than the first peak of 
the radial distribution function playing a bigger role for 
the phase behavior on the solid side 53 than on the liquid 
side. 
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FIG. 13: Test of Rosenfeld-Tarazona scaling (Eq. (J29J)) plot- 
ting potential energy versus T 3//5 on isochores for the systems 
investigated in Fig. [4] 



B. An approximate equation of state for 
generalized LJ systems 

Equations p^|) and (I2TJ|) give the potential energy and 
the virial as functions of density along an isomorph. The 
temperature variation along the isomorph can be found 
by integrating Eq. ([9]). Since 7 generally changes only 
little (Fig. [3]), in many situations it is a good approx- 
imation to assume it constant, leading to the relation 
T = p 7 T* along an isomorph, where T* is the temper- 
ature at the reference point. Taken together Eqs. (fl"9|) . 
(|2T)1) , and thus imply that from a single reference point 
(p*, T*, U*, W*) we have a prediction for (p, T, U, W) on 
the isomorph to which the reference point belongs. 

Figure 2] demonstrated that in the WU phase dia- 
gram, isochores to a good approximation are straight 
lines with slope 7, i.e., that one can write W(p, T) = 
Wo{p)+j{p)U(p, T). The only ingredient missing to gen- 
erate a full equation of state expressed as U = U(p,T) 
and W — W(p, T) is a relation between temperature and 
potential energy or virial along an isochore. 

Rosenfeld and Tarazona^ derived from density func- 
tional theory an expression for the potential energy on 
an isochore: 



U(p,T) = Uo(p)+a(p)T 3 / 5 . 



(29) 



Rosenfeld and Tarazona noted that the expression" ... 
provides a good estimate of the thermal energy (...) near 
freezing densities only for predominant repulsive interac- 
tions" , and confirmed their expression by simulations of 
IPL potentials with different exponents. Equation (|2"9")l 
has since been shown, however, to be an excellent ap- 
proximation also for several models with attraction, in- 
cluding the KABLJ mixture.— ~— This was regarded as 
a bit of a mystery, but can now be understood as a con- 
sequence of the fact that dynamics and heat capacity of 
strongly correlating liquids are well reproduced by a cor- 
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FIG. 14: Upper panel: Test of the equation of state in the 
WU phase diagram using p — 1.200 as reference density. 
Curves and lines: theoretical predictions (Eqs. (|31[) to (J36J) ) 
for three isochores, an isomorph, and an isotherm. Data 
points: simulations of the 12-6 KABLJ mixture. The param- 
eters of the equation of state were found from the two state 
points indicated by circles, giving (Wo/N,"f,Uo/N,a/N) = 
(38.4,5.17,-8.63,2.61). By comparison, the corresponding 
parameters found by fitting to the p = 1.200 isochore in 
Figs. H and [TJ are (38.3, 5.16, -8.63, 2.61), illustrating the ro- 
bustness of the fitted equation of state. Applying the fitting 
procedure to the simulated state points on the T — 0.550 
isotherm gives the parameters (38.5, 5.19, -8.63, 2.61). Lower 
panel: Same data as above in a pT diagram. 



responding IPL system^ - a consequence of the hidden 
scale invariance discussed briefly in the introduction (see 
Papers I-III). 

Figure Q~3] tests the Rosenfeld-Tarazona prediction for 
the systems investigated in Fig. 2) The prediction works 
very well, but less so for the supercooled Wahnstrom 
BL J system, for which the formation of extended "Frank- 
Kasper" clusters at low temperatures affects the temper- 
ature dependence of the potential energy^ 

Combining Eq. (1291) with our results on isomorphs and 
isochores we can construct an equation of state for gen- 
eralized LJ systems. The idea is to map any state point 
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(p, T, U, W) to the corresponding isomorphic state point 
(p*,T„, [/*, W*) on a reference isochore (where it is im- 
plecitely understood that the parameters (Wo, 7, Uo,a) 
are evaluated at the reference isochore /?*): 



T*(p,T) = Tp-~<, p = p/p* 

U*(p,T) = Uo + aiTp- 1 ) 315 , 
W*(p,T) = W + ~fU*{p,T). 



(30) 

(31) 
(32) 



From (J7*,W*) we can calculate the contribution to the 
potential energy from the two IPL terms of the potential 
(compare Eqs. (JT4J) and (|T5|) ): 



(pjT) = 3Wo-(n-37)(^ + .(T^)3/ 5 ) ; 

m — n 

C/*(p,T = — y —!—t — Uu) 

m — n 

Finally, we use isomorphic scaling (Eqs. (fTU]) and (|2U)) ) 
to go back to the (p, T, U, W) state point: 

U(p,T) = p^UUp^ + p^U^T), (35) 
W(p,T) = ^p^U^p,T) + ^ 3 u:(p,T).(S6) 



The new equation of state given by Eqs. (|33[) to ([36]) 
contains four parameters (Wo, 7, Uq, a) evaluated at a 
reference density p*. The parameters can be calcu- 
lated from U, W CyS and /3y x at a single state point. A 
more convenient way of estimating the parameters might 
be: For a collection of (p, T, U, W) state points, at least 
two of which are non-isomorphic, use isomorphic scaling 
(Eqs. (fT9t and (12010 to get the corresponding ({/*, W*) at 
the chosen reference density, and fit them to Eq. (13"2"1) to 
determine Wo and 7. Next, using the fitted value for 7, 
fit to Eq. (J2J) in order to find Uq and a. 

To put the new equation of state to a test, the pro- 
cedure described above was used to determine the four 
parameters (Wo, 7, Uq, a) for the 12-6 KABLJ system us- 
ing merely two state points. The resulting equation of 
state is shown as full lines in Fig. 1141 where it is com- 
pared to simulation results shown as data points. The 
agreement is good. The largest deviations are seen when 
changing density away from the reference density, most 
evident for the isotherm in the WU diagram and the iso- 
morph in the pT diagram. Better agreement is expected 
if the density dependence of 7 is taken into account, but 
such a correction will not be attempted here. To summa- 
rize, combining the isomorph theory with two further ap- 
proximations, a state-point independent exponent 7 and 
the Rosenfeld-Tarazona expression, a realistic equation 
of state with four parameters is arrived at. 



V. CONCLUSIONS 

Paper IV demonstrated that strongly correlating liq- 
uids have isomorphs, i.e., curves in the phase diagram 



along which structure and dynamics to a good approxi- 
mation are invariant in reduced units. The main new re- 
sults in this paper are predictions specific to generalized 
Lennard- Jones systems and the supporting numerical ev- 
idence. 

Starting from the invariance of structure along iso- 
morphs, a prediction for the shape of these in the WU 
phase diagram was presented and shown to agree well 
with simulation results. It was shown that for a given 
system, the isomorphs all have the same shape in the 
WU phase diagram; they are simply scaled versions of 
each other. Furthermore, the isomorph shape depends 
only on the exponents m and n, i.e., all systems with, 
say, 12-6 LJ interactions have same isomorphs in the WU 
phase diagram. What differs from system to system is the 
values of the density and temperature along an isomorph. 

Throughout Papers I-V the WU phase diagram was 
referred to extensively. This particular phase diagram 
has not been used much before, but it is evident that 
for strongly correlating liquids several features stand out 
when presented in the WU phase diagram. Prominent 
examples are that isochores are straight lines (with slope 
7) and the existence of the "master" isomorph discussed 
above for generalized LJ systems. 

The class of strongly correlating liquids includes, we 
believe, all van der Waals liquids and metallic liquids (but 
excludes hydrogen-bonding, covalent, as well as strongly 
ionic liquids). More simulations, as well as experiments, 
are needed to substantiate that the class of strongly cor- 
relating liquids is this large, but preliminary simulations 
of simple molecular models like2£ the asymmetric dumb- 
bell and the Lewis- Wahnstrom OTP three-site model are 
encouraging. In particular, the isomorph properties of 
more complex molecule models need to be simulated, as 
well as theoretically contemplated; the scaling properties 
of molecules with fixed bond lengths are not trivial. It 
should be noted that even a quite complex system like a 
phospholipid membrane can be strongly correlating with 
respect to its slow degrees of freedom^ the interactions 
of which are dominated by van der Waals forces. An- 
other interesting question is to study in detail the prop- 
erties of liquids where the exponent 7 varies significantly 
throughout the phase diagram, for instance the Weeks- 
Chandler- Andersen cutoff variant of the generalized bi- 
nary LJ liquid that was recently investigated om detail by 
Coslovich and Roland^ as well as Berthier and Tarjus^i 
Another system that deserves further investigation is the 
200-100 LJ system^ which we find to be an exception to 
the general observation that generalized LJ systems are 
strongly correlating - presumably as a consequence of its 
extremely narrow range of interactions. 

This series of papers investigated in depth the prop- 
erties of strongly correlating liquids. We believe to have 
demonstrated that strongly correlating liquids are sim- 
pler than liquids in general. Hopefully this insight will 
turn out to be useful in furthering the understanding of 
the properties of liquids in general. 
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